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Abstract. Starting from the concept of involution of field equations, a universal method 
04 ' is proposed for constructing consistent interactions between the fields. The method equally 

I well applies to the Lagrangian and non-Lagrangian equations and it is explicitly covariant. 

No auxiliary fields are introduced. The equations may have (or have no) gauge symmetry 



(N 



O , and/or second class constraints in Hamiltonian formalism, providing the theory admits a 

, Hamiltonian description. In every case the method identifies all the consistent interactions. 

in 

1. Introduction 

p I In this paper, a universal method is suggested for solving the problem of consistent 

i interactions in the general field theories. The method equally well applies to the Lagrangian 

or non-Lagrangian field equations, and it does not break the explicit covariance. The field 

equations may have (or have no) gauge symmetry, and/or have (or don't have) the second 

CN ■ class constraints in the corresponding Hamiltonian form in the variational case - the method 
00- 

I provides identification of all the consistent interactions, in all the instances. 
^ . The proposed method is based on the idea of involutive form of the field equations. Every 

^ ■ regular system of field equations can be equivalently reformulated in the involutive form. 

The gauge algebra of the involutive equations is more rich, in general, then the algebra 
^ ■ of the equivalent non-involutive system. In particular, the involutive system may have 
i gauge identities even if the theory does not have any gauge invariance. The consistency is 
completely controlled by the stability of the gauge algebra of involutive system with respect 
to inclusion of the interaction. 

The paper is organized as follows. In the next section, we explain the problem setting, 
and discuss previously known methods of its solution. In Section 3, we explain the notion of 
the involutive form of field equations and preview the basic idea of our method. In Section 
4, we describe the gauge algebra associated to the involutive system of field equations. 
Then, we explain the procedure of perturbative inclusion of interactions in the involutive 
systems. In Section 5, we consider the examples of interactions in massive spin 1 and 2 
models, to illustrate the general method. The Conclusion contains a brief discussion of the 
paper results. In the Appendix we derive the relation (jH]) which is used in the paper for 
counting physical degrees of freedom. 
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2. The consistency problem of interactions 

The most common view on the problem of consistent interactions can be roughly for- 
mulated in the following way. The starting point is a system of linear field equations, or 
a quadratic Lagrangian, which is supposed to be covariant with respect to certain global 
symmetry group (the most often examples are Poincare, AdS or conformal groups). The 
free field model has a certain number of physical degrees of freedom. Switching on an inter- 
action means inclusion of nonlinear covariant terms into the free equations. The interaction 
is said consistent if and only if the nonlinear theory has the same number of physical degrees 
of freedom as the original free model has had. 

If the free field equations and the interacting ones follow from the least action principle, 
the Dirac-Bergmann algorithm |1] will allow one to examine the number of physical degrees 
of freedom and to check in this way consistency of the interaction. There is an extension 
of this algorithm [2j which is applicable to general regular dynamics with not necessarily 
Lagrangian equations. The Dirac-Bergmann technique and its extensions, however, break 
the explicit covariance. Because of that, these algorithms can help to examine the consis- 
tency of particular interaction a posteriori, but they are hardly able to serve as a tool for 
covariant derivation of the consistent interactions. 

The explicitly covariant method is known for solving the interaction consistency problem 
for the gauge fields in Lagrangian theories (for introductory review see [3]). The basic idea 
of this method is that both the action and the gauge symmetry must be simultaneously 
deformed by interaction in such a way that the number of gauge transformations for the 
Lagrangian would remain the same after inclusion of the interaction as it has been in the 
free theory, though the symmetry can change. The most systematic form of this method, 
being based on the cohomological view of the problem, is known in the framework of the 
BV-BRST formalism [1], [3]. In its turn, the BRST approach to the consistency problem 
is based on the general theory of local BRST cohomology [5J. Many gauge field theories 
are known where this method either has delivered the complete solution to the consistency 
problem, or established a no-go theorem for the interactions. We mention several examples 
of such results: establishing all the consistent interactions for p-form fields [6], and the no-go 
theorem for graviton-graviton interactions [7]. There are numerous results of a similar type 
obtained in various models by various versions of this method during the two recent decades. 
Among the most recent ones we mention the work [8j, where all cubic electromagnetic 
interactions are found by this method for the higher-spin fermionic fields in Minkowski 
space, and it is proven that the minimal couplings are not admissible. 



CONSISTENT INTERACTIONS AND INVOLUTION 3 

The above mentioned method is so popular because of the exphcit covariance and alge- 
braic elegance. This method, however, is unable to provide a solution for the consistency 
interaction problem for any field theory. The matter is that gauge invariance is not the 
only cause that makes "nonphysical" some of degrees of freedom. For example, the num- 
ber of physical degrees of freedom is less than the number of fields for the second class 
constrained models, though there is no gauge symmetry. That is why, even in gauge invari- 
ant systems, the non-violation of gauge symmetry by the interaction does not necessarily 
mean the consistency of interaction. With this regard, it is relevant to mention the idea 
of introducing Stiickelberg fields to gauge the massive models, and then apply the usual 
method to introduce the gauge invariant interactions. Many works apply the Stiickelberg 
gauge symmetry idea to find the consistent interactions or examine consistency for various 
massive fields. Among the examples of this type we can mention the series of recent works 
on massive gravity [9], ^Q\, [H] (further references can be found in these articles, and also 
in [I2])- Inclusion of the Stiickelberg fields is rather an art than a science at the moment. 
No general prescription is known of doing that in such a way that could ensure the consis- 
tency of interaction just as a consequence of its consistency with the Stiickelberg symmetry. 
To examine consistency of interactions in the Stiickelberg gauged model, a non-covariant 
Dirac-Bergmann constrained analysis still remains inevitable in many cases, see for example 

In this paper we suggest to control the consistency of interaction by exploiting the in- 
volutive form of the field equations. As we explain in the next section, any field theory 
model can be equivalently formulated in the involutive form, and the result of such a re- 
formulation is termed an involutive closure. Normally, the involutive closure retains all 
of the symmetry of the original system that makes the method convenient for studying 
covariant field equations. The involutive closure of the Lagrangian equations is generally 



as Proca equations. There is no 
in the involutive closures of the 



not Lagrangian anymore, even for such a simple mode^ 
pairing between gauge identities and gauge symmetrie 
Lagrangian equations in general. In particular, even if the Lagrangian does not have gauge 
symmetry, the involutive closure of the Lagrangian equations can possess non-trivial gauge 
identities. The gauge algebra of the involutive closure of the field equations, with unrelated 



^The gauge identities are often termed as the Noether ones, because the second Noether theorem states 
the isomorphism between the identities and symmetries. For the general non-Lagrangian system of equa- 
tions, including involutive closure of Lagrangian one, there is no automatic Noether's correspondence 
between symmetries and identities. That is why, we do not use the term "Noether identity" to avoid the 
impression that it is related to any gauge symmetry. 
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gauge symmetries and identities, turns out to be more rich, in general, than the gauge 
algebra of the equivalent non-involutive system. It is the structure of the gauge algebra of 
the involutive closure of the field equations that controls the number of physical degrees of 
freedom. That is why, when the interactions are included, it is the stability of the gauge 
algebra of the involutive closure of the field equations that ensures the consistency of the 
theory. 

The idea of involution is well developed in the theory of ODE's and PDE's, and it is 
effectively applied to a broad range of the problems in various fields, see [13\. It has 
been never systematically studied, however, as a tool for controlling the consistency of 
interactions - to the best of our knowledge. In the next section we explain some of the basic 
notions related to the involutive systems which are relevant in the context of consistent 
interaction problem. Our consideration, being at the physical rigor level, skips subtleties 
of the theory of involutive systems (for a rigorous review see [I3]), because we focus at the 
issues of gauge algebra that are underdeveloped in this field at the moment. 



3. Involutive form of field equations 

By the order of a field equation we understand the maximal order of the derivatives of 
the fields in the equation. The maximal order of the equations in the system is said to be 
the order of the system. 

Definition. A system of order t is said involutive if any differential consequence of the 
order less than or equal to t is already contained in the system. 

In the theory of involutive PDE's fL3\, the definition above is normally complemented by 
more strong and intricate requirements which are inessential here, given the physical rigor 
level accepted and the problem addressed, so we adopt this simple definition. 

If the system of field equations is supplemented by all the differential consequences of 
the orders lower or equal to the order of the system, it becomes involutive. The system 
brought to involution in this way is said to be the involutive closure of the original system. 
Obviously, the original system is equivalent to its involutive closure in the sense that all 
the solutions are the same for both systems. 

If the system of field equations can be brought to the involutive closure just by inclusion 
derivatives of some of the equations, we will also consider it involutivcl, though it is not. 



In the theory of PDEs such consequences are termed trivial integrabihty conditions. These are inessen- 
tial for the count of physical degrees of freedom. 



CONSISTENT INTERACTIONS AND INVOLUTION 5 

according to the definition above. This simphfication is convenient, and it does not lead to 
any contradiction as far as the problem of consistent interactions is concerned. 

To illustrate the distinctions between the involutive and non-involutive form of field equa- 
tions from the viewpoint of variational principle and gauge algebra, consider the example. 

Example: Irreducible massive spin 1 field in d = 4 Minkowski space. The field equations 
of the model include the Klein-Gordon equations and the transversality condition for the 
vector field, 

(1) = (□-m2)A^ = 0, ord{T^) = 2 

(2) = = , ord (Tx) = 1 

The order of equations is denoted by the symbol ord. This system is involutive and it is 
obviously non-Lagrangian as it contains five equations for the four component field Afj_. 
There is no gauge symmetry, though there exists a non-trivial gauge identity between the 
equations ([1]) and ([2]): 

(3) 9^T^-(n-m2)T^ = 0. 

It is useful to introduce the generator of gauge identity L, such that the identity would read 

(4) L'^T.^O, a=(/i,±), L^ = d^, = -(□ - m') 

Define now a notion of the total order of gauge identity. Suppose the a-th component of 
the generator of the gauge identity L"^ is a differential operator of order Sa, and the order 
of the equation is t^. The order of the a-th identity component la is defined as the sum 
la = Sa + ta. The total order of the gauge identity generated by L is defined as the maximum 

(5) ord{L) = max{/a} . 

a 

In the case of the identity ([3]), the differential operator = — (□ — m^) has the order 
s± = 2, and it acts on the first order equation Ti, so the order of the transverse component 
of identity l± = 2 + 1 = 3. The order of the differential operator = 9^ is s'^ = 1, and it 
acts on the second order equation T^, so the order of this identity component is the same, 
1^ = 1 + 2 = 3. As we see, the total order of the gauge identity ([3]) is 3. 

Alternatively, the same massive vector field model is described by the Proca equations 

(6) = iKO- d^d" - 5';)A, = , ord(P^) = 2 , 
that are Lagrangian. The Proca equations follow from (II]), (E]): 

Pfi = Tfj^ — dfj_T± , 
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and vice versa, the Klein-Gordon and transversality equations follow from ([6]), 

Notice that the Proca equations, being of the second order, have the first order differential 
consequence - the transversality condition (|2]). This means that the Proca system is not 
involutive. Obviously, and (II]), (E]) are equivalent systems of equations. The involutive 
closure of Proca equations, that includes (j6]) and their first order consequence ([2]), is an 
involutive non-Lagrangian system that contains the gauge identity of the third order: 

(7) 9^P^ + m^T^ = 

One can observe that the involutive form of the spin one equations, being non-Lagrangian, 
has a gauge identity, so its gauge algebra is non-trivial. Quite opposite, the Lagrangian 
Proca equations equations are not involutive and have trivial gauge algebra, without any 
gauge identity. 

A similar conclusion as in the above example holds in general: if the Lagrangian equa- 
tions are not involutive, their involutive closure, being a non-Lagrangian system, will have 
the gauge algebra with more independent gauge identity generators than the original La- 
grangian system has had. Notice that if the Lagrangian system has a gauge symmetry, 
its involutive closure will obviously have the same symmetry, so the gauge algebra of the 
involutive closure will always have the original gauge symmetry as a subalgebra. 

Running a couple steps ahead, notice that the structure of the gauge algebra of involutive 
system unambiguously defines the number of physical degrees of freedom by the following 
formula: 

oo 

(8) ^^=J2k{h-lk-rk). 

k=0 

Here tk is a number of equations of order k in the involutive system, Ik is the number of 
gauge identities of k-th total order, and is the number of gauge symmetry generators of 
kth ordei§. In this formula, both the gauge symmetry and gauge identity generators are 
supposed irreducible. The formula is derived in the Appendix, where one can also find its 
extension to the case of reducible gauge identities and symmetries. The number of physical 

■^By the order of gauge symmetry generator we understand the highest order of the derivative of the gauge 
parameter involved in the gauge transformation of the fields. From this point of view, in Lagrangian theory, 
where the gauge identities and gauge symmetries have the same generators, the order of symmetry can be 
different from the total order of the identity. This is because the total order of identity is defined taking 
into account the order of the equations it involves, while the order of the gauge symmetry is indifferent to 
the order of equations. 
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degrees of freedom is understood here as a number of the Cauchy data needed to define a 
solution modulo gauge ambiguity. In variational case it coincides with the dimension of the 
reduced phase space. So, if the configuration space count is done, N has to be divided by 
2. 

In the example above of the involutive equations for the massive spin-1 in c? = 4, we have 
one first order equation and four second order ones, so ^2 = 4, ti = 1. There is one identity 
of the third order = 1, and no gauge symmetry. Substituting these numbers into the 
general formula ([8]) one obtains M = 1-1 + 2- 4 — 3-1 = 6 that provides the correct answer, 
as the massive vector field has 3 physical polarizations, and the physical phase space is 
6-dimensional. 

Consider one more example illustrating relation (|8]): a scalar field x) in two-dimensional 
space subject to a pair of the second order equations 

(9) Tt = d^^{t, x) = 0, = dl^{t, x) = 0, ord{Tt) = ord{T^) = 2 . 

There are no differential consequences of the second or lower orders, so the system is in 
involution. There is the fourth order gauge identity, 

(10) d^T, - dlTt = . 

With two second order equations {t2 = 2) and one forth order identity (^4 = 1), relation ([8]) 
brings zero as the number of physical degrees of freedom. Let us directly check that it is a 
correct count. The general solution reads 

(11) (j){t,x) = Axt + Bx + Ct + D, 

with A,B,C,D being arbitrary integration constants. No arbitrary function of a; or t is 
involved in the general solution. This means that the system, being 2d field theory, has no 
local physical degrees of freedom. 

Let us comment on formula (IHl). At first, we notice that the relation is valid for involutive 
equations. So, prior to applying this formula to a non-involutive system, one has to take 
an involutive closure. The involutive closure can be taken in an explicitly covariant way 
for covariant equations, it does not require any (3 + l)-splitting, unlike the Dirac-Bergman 
algorithm. Also, it is important that the involutive closure can be found for any system, 
be it Lagrangian or not, and it does not require any special (e.g., first order) formulation. 
The orders of equations, symmetries and identities can be easily found, so formula ([8]) 
provides a simple tool to covariantly control the number of physical degrees of freedom. 
The second peculiarity of this formula is that it is somewhat counterintuitive: it involves 
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neither the number of fields in the system, nor is it sensitive to the number of zero order 
gauge symmetries, identities and equations. 

The problem of computing the number of the physical degrees of freedom has been thor- 
oughly studied in the theory of involutive systems pL3j, though the answer has been known 
in a completely different termqj, not explicitly appealing to the total orders of equations, 
identities and symmetries. In the context of the problem of consistent interactions, where 
the structure of the gauge algebra is the principal object to study, it is important to control 
the degrees of freedom in terms of the gauge algebra constituents. In the Appendix we 
deduce the formula (IE]) proceeding from the definition accepted in the theory of involutive 
systems and based on the concept of strength of a system of equations. This concept was 
introduced by Einstein when he counted degrees of freedom in General Relativity [H] and 
it has been further developed in many works (among which we mention [15], |16].|17].|18]) 
and related to the count by Cauchy data [13]. 

Notice that for Lagrangian equations, whenever they are involutive from the outset (in 
this case, the gauge identity generators coincide with the gauge symmetry ones), the receipt 
dH]) for the degree of freedom count takes a special form fl82|) . which has been well known 
before [19]. Notice another special form of field equations where the asymmetry may occur 
between gauge identities and symmetries: the unfolded formulation of the higher spin fields 
(for review see [20]). This method utilizes the involutive form of the unfolded equations, 
and it also benefits from the fact that all the equations, symmetry and identity generators 
are of the first order. The unfolded formalism involves, however, infinite number of field 
equations, symmetries and identities. Formula ([8]) can not be immediately applied to the 
unfolded systems because all the numbers in ([H]) are supposed to be finite. In this case, the 
method of (T_-cohomology [21], [22], [23] provides a tool for the degree of freedom count. 
The (T_-cohomology method allows one to pick out a finite involutive subsystem such that 
the unfolded system will be its infinite jet prolongation. The degree of freedom count in 
the finite subsystem, being made by the formula ([8]), delivers the answer for the number of 
degrees of freedom in the complete unfolded theory. 

Let us formulate now the key stages of the procedure we suggest for constructing consis- 
tent interactions, given the original free field equations: 

(1) The free system is to be brought to the involutive form. 



^The analysis is made by means of the theory of Cartan's differential systems, and the answers are 
formulated in terms of the Hilbert polynomials |13j . In principle, this way of counting physical degrees 
of freedom is sufficient, though it seems inconvenient in the context of relativistic field theories because it 
requires cumbersome and not always explicitly covariant computations. 
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(2) All the gauge symmetries and identities are to be identified in the free involutive 
system. 

(3) The interaction vertices are perturbatively included to comply with the three basic 
requirements in every order of coupling constant: 

(a) The field equations have to remain involutive; 

(b) The gauge algebra of the involutive system can be deformed, though the number 
of gauge symmetry and gauge identity generators remains the same as it has 
been in the free theory; 

(c) The number of physical degrees of freedom, being established by relation (|8]), 
cannot change, though some of the involved orders can. 

This procedure ensures finding all the consistent interaction vertices, for any regular system 
of field equations. 

4. Gauge algebra of involutive systems 

As it has been already explained, if the Lagrangian system of field equations is not 
involutive, its involutive closure will be non-Lagrangian. It is the structure of the gauge 
algebra of the involutive closure, not the original system of equations, that controls number 
of physical degrees of freedom. That is why, one has to study the gauge algebra of dynamics 
in its non-Lagrangian involutive form even if the Lagrangian exists. The exception is the 
case where the Lagrangian equations are involutive from the outset. In this special case, 
the known methods [3j work well, being insufficient for general Lagrangians. Notice once 
again that the non-involutive Lagrangian equations are not exceptional - many models in 
physics are of this class, e.g. massive fields with spin. This leads us to consider first the 
gauge algebra of the non-Lagrangian dynamics. 

The gauge algebra of the general (not necessarily Lagrangian) system is known in the 
same details as in the Lagrangian case, and the corresponding BRST complex is also well 
studied [21], |25] that allows one to systematically control all the compatibility conditions. 
Below, we provide a simplified description of the gauge algebra, without resorting to the 
corresponding cohomological tools and leaving aside the higher compatibility conditions, as 
these are less important in the context of interaction problem. 

4.1. Algebra of gauge symmetries and identities in general field theory. It is 

common to consider general structures of gauge algebra by making use of the condensed 
notation, and we will follow this tradition as it is convenient for presenting the general idea 
of the method. In this notation, the fields are collectively denoted by with i being the 
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condensed index that includes all the discrete indices, and also the space-time points. For 
example, the vector field A'^i^x) is indexed hj i = {fJ.,x). Summation over the condensed 
index implies integration over x. 

In the condensed notation, any system of field equations reads 

(12) 7;(0) = O, ord{Ta)=ta, 

where a is a condensed index, and is understood as a function of the fields and their 
space-time derivatives up to some finite order ta. The discrete part of the condensed index a 
labeling the equations is different, in general, from that of the condensed index i numbering 
the fields. For the Lagrangian equations, i coincides with a, though this is not true if the 
involutive closure is considered instead of the original equations. For example, the involutive 
closure of the Proca equations includes both the original Lagrangian equations ([6]) and the 
transversality condition (|2]), so the indices belong to the different sets: i = (/i, x), and 
a = (/i, For the regular field equations, the order ta depends only on the discrete 

part of the index a, not on the space-time point. 

The general field equations can enjoy gauge symmetry transformations 

(13) 5,0' = e"K(0), 5eT,(0) 1^^0 = 0, Ve" , ord{R^) = < oo , 

where the gauge parameters e"" and generators Ra{(f)) are understood in the sense of con- 
densed notation, i.e. the summation over a implies integration over x. For example, in 
electrodynamics, 6eA^{x) = d^e{x), and hence i = (/i,a;), a = y, R^{x^y) = d^5{x — y), so 
that 5eA^{x) = f dyR^{x,y)e{y). Locality of the gauge symmetry implies that the gauge 
generators R^ are the differential operators of finite order, denoted r^, with coefficients 
depending on the fields and their derivatives. 

The condition (fT2|) defines the on-shell invariance of the equations that off-shell reads 

(14) RMdM^) = KM)U4>) , 

where the derivative di is understood as variational and the structure coefficients U^ai^) 
are supposed to be regular on shell. 

The gauge identities can also take place for the general field equations, being not neces- 
sarily related to the gauge symmetries 

(15) L^(0)r,(0) = 0, ord{LA) = la. 

The gauge identity generators La are supposed to be local differential operators. The total 
order of the identity ord{LA) is defined by the order of the differential operator and the 
order of the equation it acts on as explained in Section 3 below relation (jl]). 
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The gauge symmetry and gauge identity generators are considered as trivial whenever 
they vanish on shell, that is 

where p and C, can be arbitrary local differential operators of finite order with the coefficients 
depending on the fields and their derivatives. 

The sets {L"}^ of ^^e gauge symmetry and gauge identity generators are supposed 

to be complete. The completeness means that any other generator of gauge symmetry or 
identity, satisfying or (ITSl) . must be a linear combination of the generators from the 
given set modulo the trivial ones ( fT6l) . 

Let us discuss now the equivalence relations for the systems of field equations. Two 
systems of field equations are considered as equivalent if they are related by a locally 
invertible change of fields and/or by locally invertible linear combination of the left hand 
sides of the equations. The admissible class of changes of fields reads 

(17) 0^^0'^ = 0'^(0,50,920,...), 

where the existence is implied for the inverse change belonging to the same class, i.e. the 
original fields have to be unambiguously determined by ( fT71) as functions of the fields (p' 
and their derivatives up to some finite order. 

For example, consider the system of vector and scalar field. The change A'^ = — 
d^(f), 0' = is admissible as the local inverse change exists. 

Given a set of fields, the admissible class of equivalence transformations for the field 
equations reduces to the linear combining with invertible coefficients: 

where the elements of the transformation matrices K and are the differential operators 
of finite order. 

Assuming the completeness of the generators of gauge sjnnmetries and gauge identities, 
one can derive the following consequences from the relations ( IT^ . (fTSjl : 

(19) Rid.ir^ - I^^d.K = U2f,Rl, + Wl^lTa ; 

(20) IV,d,L\ = U^^Ll + W^\T, , W:\ = -Wl\ , 

where ?7, W are some structure functions. These relations have further compatibility condi- 
tions involving higher structure functions (see for details [23], [25]). The existence of all the 
higher structure functions and their locality have been proven in [26] under the condition 
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that the generators L, R and the structure function U involved in (HM are all local. The 
corresponding existence theorem for Lagrangian theories has been known long before [27] . 
So, with the existence theorem, to ensure consistency of the field theory ([T2|) . it is sufficient 
to fulfill relations (JT4l) . ( IT5l) with some differential operators R, Z, U of finite order. 

4.2. Gauge algebra and perturbative inclusion of interactions. Consider involutive 
system of free field equations 

(21) tW(0) = O, orrf(Ti°)) = . 

As the free field equations are supposed to be linear, the generators of gauge symmetries 
and gauge identities are the differential operators with field-independent coefficients. With 
this regard, relations (JT4l) . (|T5l) in the free theory should have identically vanishing on-shell 
terms: 

(22) i?(°)^a.Ti°)(0) ^ , orrf(i?i°)) = ri°) ; 

(23) 4°)"T,(0)^O, ord{Lf) = l^^\ 

Given the orders of the equations, gauge symmetries and identities, the number of physical 
degrees of freedom in the free model, J\f^'^\ is defined by ([8]). 

Perturbative inclusion of interaction is understood as a deformation of the equations, 
identities and gauge symmetries by nonlinear terms. 



(24) Ti°) -> = + ^Ti^) + /Ti^) + . . . , 

(25) R'^^' ^ K = R^a^' + 9Rl^^' + + • • • , 

(26) Lf'' ^ L\ = Lf" + gL^l^"" + g^L^l^"" + .... 



Here is a coupling constant considered as formal deformation parameter, generators L\ 
and -Ra^* are linear in fields and their derivatives; Ti^'', L^''", and i?!^^* are bi-linear, etc. 
Notice that in each order of the deformation, the orders of equations, identities and sym- 
metries can never decrease. 

Now, we can give a more specific formulation of the consistency conditions for the in- 
teractions than the general explanation in the end of Section 3. The consistency of the 
interaction is provided if the three conditions are fulfilled: (a) the system remains involu- 
tive at each order in g\ (b) the deformations do not break the gauge algebra generated by 
relations ([HD, (fT5|) . though the structure functions can change in as well as the higher 
relations; (c) the orders of the equations, symmetries, and identities may increase, though 
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the overall balance established by relation ([8]) cannot change, i.e., it is required J\f^^^ = J\f 
in every order in g. The conditions (a) and (b) provide algebraic consistency of the system 
with perturbatively included interactions, and (c) ensures that the interacting system has 
the same number of physical degrees of freedom. 

Let us elaborate on the perturbative procedure of the interaction inclusion. Suppose 
we have taken the most general ansatz for Ta^\ ord{Tjf' + gT^^) = tii \ that does not 
break involutivity, so (a) is fulfilled. Substituting this ansatz into relations f lT4|) . f lTSj) and 
considering that in the first order in g, we find the following relations between T^^\ R^^\ L^^^: 



(27) R^^^%T('^ = Ui'^'Tt^ - i?«^5.Ti°) , 

(28) l'^Tt^'^ + L^A^Tj^""^ = 0. 

These relations impose nontrivial restrictions on the first order interaction. 

The first relation means that the free theory gauge transformation has to leave the 
first order interaction on-shell invariant modulo linear combination of variations of free 
equations. "On-shell" hereafter means on the free equations. 

The second relation means that the free gauge identity generators must leave the first 
order interaction on-shell invariant. Notice that even if the model has no gauge symmetry, 
the involutive closure of its equations can have non-trivial gauge identities. This means that 
the conditions (|28ll are essential for consistency of interactions even in the systems without 
any gauge invariance. If the equations are Lagrangian and involutive from the outsej^, 
relations fl28|) are reduced to the the on-shell gauge invariance of the cubic vertices in the 
Lagrangian. For the involutive Lagrangian equations, the generators R and L coincide, and 
the relations ( 127|) follow from ( l28l) in this case. For general system, including the involutive 
closure of the Lagrangian equations, the relations fl28|) are not necessarily connected with 
(1271) . The relations ( 1271) are always first examined in the Lagrangian case (see [3]- [8]) to 
check the first order consistency of the interaction in Lagrangian dynamics. As is seen 
from the explanations above, if the Lagrangian field equations are not involutive, the first 
order consistency of interaction requires to independently impose the extra conditions ( l28l) 
on the vertices, because these are not necessarily connected to the gauge symmetry of the 
Lagrangian. 

Given a free model, the solution does not necessarily exist in any theory for the first 
order consistency conditions ( 127|) . ( 128|) imposed on the first order interactions T^^^ and the 



^As it has been already noticed in Section 3, there are many non-involutive Lagrangian equations being 
of a considerable interest in physics. 
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corresponding first order corrections to the gauge identity and gauge symmetry generators 
L^'^\ R^^\ If a solution exists, it can be explicitly found as the system is linear. 

The solutions for interactions are considered modulo ambiguities related to the equiva- 
lence relations (JT8i) . ( IT71) . In particular, a nonlinear change of fields in the free equations is 
not considered as an interaction as well as a linear combination of the free equations with 
field-dependent coefficients. 

If the order t^a^ increases because of the first order interactions, the orders of gauge 
identity and symmetry generators should also increase in a corresponding way to have the 
same number of physical degrees of freedom ([H]). If a solution to ( !27j) . ( l28l) exists with the 
correct A/", one can proceed to the next order. 

In the second order in g, the basic relations of the gauge algebra (fT^ . (fT5l) read 

(29) i?(°)^9,Ti^) + i^W^a.Ti^) + i?i^)^9,Ti°) = t/(?T(^) + [/1?t(°) , 

(30) Lf^Ti^^ + L^^Ti^) + 4')'^Ti°) = . 

In the first instance, these relations represent further compatibility conditions for the first 
order interaction. Let us explain that in the case of relation fl30l) . On substituting into 
( |30ll the expressions for L^^^'^,Ta^^ previously derived from ( |271) . (1281) . one has to get a 
combination of the free theory gauge identity generators modulo free equations. This 
requirement is not automatically fulfilled for any interaction derived from fl271) . fl28|) . In 
some models it can be even possible that these relations are inconsistent. In this case, 
one arrives at a no-go theorem for the interaction. So, the second order relations (!29|) . 
fl30|) provide an additional selection mechanism for the first order interactions. If this 
filter is passed by the first order interactions, then relations (1291) . (1301) can be viewed as a 
consistent algebraic system of linear equations defining the second order contributions to the 
equations, gauge symmetry and gauge identity generators: Ta^\ Ra^\ l!^"^"" . The solution for 
the second order interaction is to be considered modulo the equivalence relations (IT71) . (|T8l) . 
In particular, the nonlinear changes of fields or combinations of the lower order equations 

(2) 

with field dependent coefficients are not considered as interactions. If the solution for Ta 
involves the field derivatives of a higher order than T^^ and then the orders of the 
gauge symmetries and gauge identities have to increase in the corresponding way to provide 
the same number of physical degrees of freedom according to relation ([8]). 

On substituting the second order interactions into the third order expansion terms of the 
relations (fl^ . (fT5|) one arrives at the relations that represent the consistency conditions 
for Ta'^\ Ra^\ L^^^'^ . This is much like the relations ( l29l) . ( 130|) work for the previous order 
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equations and generators. Again, there can be either inconsistency found at this stage, or 
one derives the third order interaction, and the procedure repeats in the next order. Three 
different scenarios are possible for further development of the iterative constructing the 
interactions. The first is that the iterative analysis of the expansion in g of the conditions 
(fT^ . (fT5|) will terminate at certain order because of contradiction. This results in a no 
go theorem. The second is that starting from certain order all the interactions become 
trivial. This results in a polynomial interaction. The third option is that the procedure 
results in nontrivial consistent interactions in every order. This leads to a non-polynomial 
interaction. 

In the case of involutive Lagrangian equations, this procedure reduces to the commonly 
known method of inclusion interactions between gauge fields (see [3] for a review). It has 
been already mentioned that the involutive closure of non-involutive Lagrangian equations 
is not Lagrangian anymore. Because of that, the Lagrangian method does not apply to this 
case, whereas the method of this section still works well, as well as for any other involutive 
system of field equations. Our method exploits the same general idea as the Lagrangian 
gauge approach: to include interactions by a consistent deformation of the equations and 
the related gauge algebra. The main distinctions are related to the fact that the general 
gauge algebra of involutive system involves gauge identities ( |T5|) independently from gauge 
symmetries, and the involutive form of equations allows one to effectively control the number 
of physical degrees of freedom. 

5. Examples: Consistent self-interactions of massive fields of spin 1 and 2 

In this section we illustrate the general procedure of perturbative inclusion of consistent 
interactions described in Section 4 by the examples of self-interactions for massive fields of 
spin one and two. Applying this method we find all the consistent interaction contributions 
(without higher derivatives) to the field equations of the second order in fields. For the 
corresponding Lagrangians this would correspond to the cubic vertices, though we find 
that some of the admissible interactions do not follow from any Lagrangian. 

5.1. The massive spin 1 in d=4. As it has been explained in previous section, there are 
equivalence relations for the involutive field equations, so one can choose various represen- 
tatives from the equivalence class of the free equations. For the spin 1, this choice is not 
unique either, as it has been explained in Section 3. We choose the Proca equations and 
the transversality condition as free involutive equations for the spin 1, 

(31) T^°) = d^'F,, - m^A,, Tf^ = A, , ord(TW) = 2 , ardirf^) = 1 . 
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In this section, we adopt the following agreement for the strength tensor and its dual: 
F^y = d^Ay — di^A^, Ffj_y = \e^yapF'^^ . The choice of free equations in the form ( l3Ti) 
is slightly more convenient than the other equivalent options, e.g. ([1]), ([2]), because the 
second order equations in the system (j3Ti) are Lagrangian that makes it easier to check the 
consistency of the interactions. 

Equations fl3Tl) admit the gauge identity whose generator reads 

(32) L(°)^ = L(o)^ = , L^o^^T^ + L^^rf = , ord{L^^^) = 3. 

So, the involutive form of the massive spin-1 free field equations includes four second order 
equations and one of the first order together with the third-order gauge identity between 
them. (The general definition is provided by relation ([5]) for the total order of the gauge 
identity.) 

The next step according to the general procedure of perturbative inclusion of interactions, 
as described in Section 4, is to take the most general covariant ansatz for the first order 
correction to the field equations that does not break involutivity. Let us assume that no 
higher order derivatives of the fields are included in the interactionqj. Then the most general 
ansatz reads 



Tj^^ = A^ [p^F,^ + psi^M") + P2 {A^d^A^ - t|°^ A^) + O^A^ [p,dpF^,+ 

+P^dpF^^ + p^dJn^^ + {pjdo^Apd^A^ + p^dpA^d'^AP + peF^^jF"^) + 
+CiA^Tf + C^d^A^d^Tf^ + C.d'^A^d^Tf^ + UF^d^Tf^ + Cs^M (^f )' + 
+C&T^a^d^A^ + CiT^^'^d^A^ + CsF%T^a^ + C^Tf^T^^^ , 

t|') = pioF^^F^p + pnd^A^d^Ap + puF'^^F^p + pnm^A^A^+ 

+ (Cio + 4)(^r^' 



(33) 



This assumption does not restrict generality. One can see that the inclusion of higher derivatives would 
inevitably increase the number of physical degrees of freedom, as it is defined by relation (|8]). We omit 
the detailed proof of this fact, though the simple evidence of that can be easily seen. If, for example, the 
third order derivatives are included into tI^^\ there will be four equations of the third order, so the positive 
contribution to Af will increase by 4. There is only one gauge identity, so its total order should increase at 
least by four to compensate that. To achieve such a growth of the order of the identity, one has to raise 
the order of Ti. As the order of the scalar equation raised, this will again increase Af with no way to 
compensate the latter growth of the order. 
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The vertices with ten ^-coefficients are trivial as they are reduced to a hnear combination 
of the free equations with field-dependent coefficients. Inclusion/exclusion of these vertices 
gives an example of the equivalence transformation (1181) with 



(34) 

(Kc);t = CiA, + [C2d,A" + C-sd^A, + CF%) + 2C,d,Tf\ (Kj^ = . 

The inverse transformation has the form (ii"^)"^ = K^(^ + 0{(). We keep these terms to 
simplify the control of trivial vertices in the next order of interactions. Besides the trivial 
terms, the most general covariant quadratic ansatz (133|1 includes a 13-parametrer set of the 
non-trivial contributions with the coupling constants p. 

Substituting the ansatz (l33ll into the structure relations (1281) . we obtain the following 
consistency conditions for tI}^\t^^: 



(35) 



= (f + m\pr + Pio))F"^Fa/3 + (f + m\2p, + Pu))f''^F^^+ 
+ (p2 + m\2ps + 2p7 - p5 + pn)^d^A''daAf^ + (^pi + m^pn)Ac,A°+ 
+ + Spejc^^F^^^^F"^ + 2p7d^d^A^dWA^ + 2psd^d^A^dWA^+ 

+p,A'^Ti'^ + d"Af^(2prd^Tp + {2p, - p,)dpT!^^ + 2{p, + p,)d^dpTf^) + 
+ApeF'^^daTP + 0(C) = (mod T(o)) . 

The C-terms vanish on-shell, and for this reason, the parameters ( remain arbitrary at this 
stage. The consistency requirement ( l35l) imposes seven conditions on thirteen interaction 
parameters p: 

P7 = 0, ps = 0, p9 = 4p6, pio = -7r^5 Pii = P5 ^ , 

_ P3 o _ P^ 
Pl2 — —7^ — o ~ ^P65 Pl3 — 2 ■ 



Obviously, six parameters pi , . . . , pg remain arbitrary, while the seven others are fixed by 
these relations. Having the consistency conditions (135|1 fulfilled, we arrive at the following 
six-parametrer set of vertices: 
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(37) ; 

T^' = -5 (^^"^^«. + -^^«^") + (P5 - ^) d^A-d.A,- 
-(£,^-2pe)F-^.. + 5(T(0))^ (modC). 
The corresponding contributions to the gauge identity generators read 

(38) L^^> = -piA" + (^ps^M" - 4p6F"^) = . 

Consider now the problem of compatibihty of the first-order interactions at the next order. 
Following the general prescription of Section 4, we have to substitute the first-order gauge 
identity generators and equations L^^^T*^^) obtained above, into relations fl30l) and examine 
their consistency. We have 

^ (l^-a^A^ - A^TfA - p.A-d'^A^ (pAd^hu + P,d^F^,+ 

(39) ^ , ^2 ^ \ ^ ^ , 
+4p6a„F^, j - pip^A-d,{F^^F^p) + [p^d-A^ - Ap^F^'^) d^T^;'^ 

As is seen from (!30|) the first order interactions L^^\T^^\ having the form (139!) with six 
parameters involved, will be compatible in the second order if there exist functions t!^'^ , T^^ 

(2) (2) 

such that ord(Ta ) < 2, ord{T\_ ) < 1 and the following conditions are fulfilled: 

(40) 2,(i)°T^i) + 5"Tf ) + ^2yj2) = (mod T(°)). 

On substituting fl39l) into fj40|) . one can find that no obstructions occur to the existence of 
Ta \ with appropriate orders of field derivatives. For example, we can always take 

^i'^ = - {p5d^A^ - 4p6F/) Tj^^ + piA^9"A^(p4F,^ + P5F„^) + 

(41) +P1P6 {AA^d^A^F^f, + A^F^^hp) , 

Tf' = ^ {^A-d^A^ - A^Tf)) - ^^d.A^d-A^ [p,h. + P,F^.) • 

This means that the six-parametrer set of the first-order interactions ( 1371) . being the general 
solution to the first-order condition 0281) . admits a consistent extension to the second order 
without any restriction on the parameters p. We will not further elaborate here on the most 
general interactions of the higher orders, although the method allows one to study the issue 
in its full generality in any order, as it can be seen from the first-order example. Instead, we 
will just notice some special cases, where the perturbative procedure of interaction inclusion 
can be consistently interrupted already at the second order level. 
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(2) 

At first, notice tliat if tlie parameters p are cliosen in such a way that = and 

(2) 

= 0, the identity will be consistent without higher order contributions, i.e., with 
n > 2. The corresponding first-order vertices are called self-consistent. 
Two special combinations of the parameters are possible that define inequivalent self- 
consistent first-order interactions: 

(1) pi , P2 , Ps 7^ 0, P4 = Ps = Pe = that corresponds to the equations at most 
cubic in fields. A further specialized option pip2 = leads to the at most quadratic 
in fields interactions; 

(2) P2 , Ps , P4 7^ 0, Pi = Ps = Pe = results in quadratic interaction. 
In the first case, the corresponding equations read 

= d-F^^ - m^A^ + A» (^p,F^^ + p^F^^^ + p^ [a'^O^A^ - Tf^A^ , 

= d'^A^ - ^ {^-F-^F^, + m'A^A-) - [d^A^d^Af, - {T^^y 



(42) 



^F-^hp + ^ { -A-d^A^ - A'Tf^ 



2m? 

L"" = d^ - p^A'^, = rr?. 

The second item results in a different self-consistent interaction of the first order 
= d^F^^, - m^A^ + P3A»F^^ + p^ [a^O^A^ - Tf A^) + p^d^A^d^h^.. 

(43) = d-A^ - ^ [d^A-d^A, - (tW)' ) - i^.F'^'h,. 

= 9", = m2 . 

These two different quadratic interactions, being self-consistent as such, can be comple- 
mented by the higher order interactions. The more general quadratic interactions ( l33l) 
need cubic corrections to ensue consistency. Though such corrections exist, as we have 
explained above, they can be inconsistent in the next order of interaction. 

Notice that the three-parameter sets of self-consistent interactions ( H2l) . ( H3ll . being the 
most general in this class, are not necessarily compatible with variational principle, though 
the free theory admits Lagrangian formulation. One can see that only the one-parameter 
family of the vertices ( H2|) . fH3|) is variational. It is the case of p2 = —pi = g and the other 
p's and C's vanishing. The corresponding Lagrangian reads 

C = + £(1) = -^F^^F^^ - ^A'^A^ - ^d'^A^A^Ap . 

All other vertices of f l42|) . fH3|) do not follow from variational principle. This demonstrates 
that the general class of the consistent interactions can be much broader than that of La- 
grangian ones. This fact can also mean that some of the no-go theorems for the interactions. 
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known in the Lagrangian framework, may be bypassed if the requirement is relaxed for the 
vertices to be variationaL 

5.2. The massive spin 2 in d=4. The irreducible spin-2 massive field theory can be 
described by a traceless, symmetric, rank-2 tensor field /i^,^ subject to the Klein-Gordon 
equations and the transversality condition 



(44) rW^(n-m2)V = 0, rW=9^V = 0, ord{T^'J) = 2, orrf(TW) = l. 



These equations are involutive as there are no low order differential consequences. Unlike 
spin 1, the equations are inequivalent to any Lagrangian system formulated in terms of 
the original irreducible field. The Fierz-Pauli Lagrangian [28] that involves auxiliary scalar 
field and traceless tensor hf^i, leads to the equations that are equivalent to (144|) . The Fierz- 
Pauli equations (FPE) are not involutive. Their involutive closure is not Lagrangian and 
it has a more complex structure than the equations (jHj) formulated without any auxiliary 
field. A similar picture is observed for all the higher-spin massive fields. The Lagrangian 
formulation due to Singh and Hagen [29] needs auxiliary fields, that makes the system non- 
involutive. The involutive closure of the Sing-Hagen equations is not Lagrangian anymore, 
and is more complex than the system of Klein-Gordon equations and the transversality 
condition for the traceless tensors. The aim of this subsection is to demonstrate by the 
example of the spin-2 field that the minimal formulation of the irreducible field equations, 
involving just the mass shell and transversality conditions, is sufficient for iterative con- 
struction of consistent interactions. Though this formulation is not Lagrangian, it admits 
quantization and can enjoy all the other advantages of Lagrangian formalism, including 
Noether's correspondence between symmetries and conserved currents. The matter is that 
the model admits a Lagrange anchor. As is known, the Lagrange anchor p5], being identi- 
fied for not necessarily Lagrangian field equations, allows one to path-integral quantize the 
theory [25], ^U\, [31], and also to connect symmetries with conservation laws [26], [32] . 

Prior to seeking for consistent interactions, we have to identify the gauge identity and 
gauge symmetry generators for the free field equations f l44|) . The model has no gauge 
symmetry and there exists four third-order gauge identities. The generators are given by 



(45) 





(46) 



lW^^tW + L^o)u^t!;') ^ , ord(L(''))„ = 3 
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Following the general procedure of Section 4, to switch on the first order interactions, 
one has to find the quadratic vertices T^iJ^^ri^^ such that the identities ( 128|) hold with 
ri?,Ti°\LL°)'^^LL°)^ defined by dH, gS]). 

We do not study the most general case, restricting the quadratic vertices by the ansatz 
with at most two derivatives in every termj^ 

where we used the following abbreviations: 

Notice that the trivial (on-shell vanishing) terms are omitted in (H7j) . All the vertices are 
identically traceless. Following the general procedure of Section 4, we substitute the ansatz 
f j47l) into the relations (l28|) and examine consistency. 



(48) -(P9 + Pi)d»hf'^d^TP + p6/^°^9,9^ri°) - p,dPh\d^Tf^ + pjh<-Pd,T^^)) + 

+p,h'^^d,d^Tp + p.h'^^d^T^'J - (p9 + p,)m'h^X'^- 



where 

Qt,u = d^h^^dahpy (^p9 + p3 + pi- 2p2^ + d^h°'^dahf}u (p4 + Pe) + 
(49) +^r]^, {dhf - P3 + 2p7 - 8pio) + ^r/^. ' (2pi - P4) + 

+w?h'^^hau {p9 + P5 + Pi - P2) + ^ni^Vf^uh^ {^Ps - P5 + P7 - 4p 



In the most general case, the terms can appear with two second-order derivatives and four derivatives 
in totaL The interactions of this type are usuaUy considered abnormal, and we do not study them here. 
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At first order, the consistency requires the expression (148|) to vanish modulo T^^\ Each 
term in (B9|) being independent (modulo a conserved tensor), so we arrive at the system of 
equations restricting the interaction parameters 



(50) 



P9 + P3+ Pi- 2p2 = 0, Ps- P3 + "^Pl - 8pio =0, P4 + P6 = , 



P9 + P5 + Pi - P2 = 0, 2ps - P5 + P7 - 4pio = 0, 

The solution to this system reads 



2pi - P4 = . 



(51) 



7 

PS = -3, 

P5 = P2 - 7 , 



P9=1-Pl, P3=2p2-7, Pi = '2p 



I5 



7 P2 
P6 = -2pi, p7 = 4po--, Pio = Po-^, 



where 7, po, Pi, P2 are arbitrary constants. The term with the coefficient p2 can be considered 
as trivial, because it is generated by the diffeomorphism in the space of fields: 



Finally, we conclude that the set of non-trivial vertices for the massive spin-2 equations 
may involve at most 3 parameters. The consistent quadratic vertices are given by 



m 



+ Pi 



(0) 



(52) 



+ 



+P2 
+P0 



+ 



T^^ = d^^ (p2/i„^/i^ + (po 



P2 

4 



Notice that two parameters 7 and pi are associated with the conserved currents that do 
not contribute to the transversality condition 



(1) 



(53) 



The first-order deformation of the Noether identity generators is given by 
^''('^^M^ ■ ) + (pi - l)d^''h^\ + Ap^d^h^P + pi^f.'TW'^) , 



r(l)°/3 



9^ £ 
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where the round brackets mean symmetrization in corresponding indices and the following 
notation is used: 

/ 1 ^ 

i^tu = Pi (5.^% + 'inh^dp - 2d'^h,, - h^d, + -r^^,h^-dpj + 

+l{d,h\-d,h-^y 

As it follows from the general requirement fl3UI) . the next order consistency of the vertices 
fl52|) is only possible under the following condition: 

(55) ^'^"'^^ii^ + L^a^'T^'^ = d^Qc, - m^Tf (mod T^) , 

where 



(56) 



^afi — 'Jfi-'-a -'-afi ■ 



Relation (H9|) is the compatibility condition for the second-order vertices. We will seek for 
the solutions to these equations that obey conditions 



ord{Q 



ordiTf^) = 1 . 



These are obviously consistent with the number of physical degrees of freedom. The ad- 
missible choice is 



n — rOlPrp(l) rp{l) I 



-r -r yPi - 1) [h^ir;,'^ - d.T^) + d.h^j'^^ 

To get Q we integrate by parts and take into account the identity fHSl) : 

The on-shell vanishing gauge identity generators and di^T^^^^^ can not affect the relation 
([55]) . so the question left reads: is it possible to represent the remaining term /i"^5,T^J^^ m 
the form fl55l) for some Q? If the answer is affirmative, one should try to express Q in the 
form (1561) with on-shell traceless and symmetric tensor T^J . If the answer is negative, an 
obstruction for the second order vertex appears. We have 



(57) 



Apod^h'^^Ti'^ = d>^ ApoUpo - Ddyh^^dpKr + ^PoPiduh^h^'^drK^ 



-SpopM^^h^^dMa^ + . 



mo 



d TW) , 
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where the dots denote the terms of order 1 that may be included into . There is a term 
of order 2 in the r.h.s. of ( 1571) 

that does not reduce on shell to a total divergence. It cannot also be absorbed by deforma- 

(2) 

tion of the transversality condition T/i because of the order restriction. This means that 
the vertices are inconsistent in the class with ord{Tjf') = 1 unless poPi = 0. Therefore there 
can be at most two 2-parameter families of consistent interactions in the considered class. 
This seems matching well the fact that the massive gravity admits a 2-parameter family of 
Lagrangians [33] that are consistent from the viewpoint of Hamiltonian constrained analy- 
sis |3l]. The detailed comparison, however, with complete nonlinear equations of massive 
gravity [33] is not straightforward as the vertices for the field equations (jHj) are deduced 
for the traceless tensor /i^,^ without any auxiliary field involved. The massive gravity equa- 
tions involve the tracefull tensor. In the free limit the massive gravity reduces to the FP 
equations. In the FP theory, the trace vanishes on shell, leading to the equations (jHj) 
for traceless tensor. For the nonlinear equations of massive gravity, the explicit on-shell 
exclusion of the auxiliary field is unknown at the moment, though one should expect it is 
still possible (if it was impossible, the theory would have a different number of degrees of 
freedom). Unless the auxiliary field remains involved in the equations, the corresponding 
vertices can not be immediately compared with the ones in the equations formulated with- 
out auxiliary fields from outset. The goal of this section, however, is not to perturbatively 
re-derive the massive gravity vertices, but to exemplify the involutive technique of finding 
interactions without a direct use of Lagrangian, making use of a minimal set of fields. 

6. Concluding remarks 

Let us briefly discuss the results. In this paper, we propose a new method that allows 
one to examine the consistency of interactions for the general field theory models, be they 
Lagrangian or non-Lagrangian. The method also provides a technique for perturbative 
identification of all the admissible interactions, given a free field model. The method re- 
quires first to bring the field equations to the involutive form. Notice that the involutive 
closure is always non- variational for variational non-involutive equations. As far as the field 
equations are brought to the involutive form, the gauge algebra is to be identified for the 
equations. The consistency of the gauge algebra is examined by tools similar to those based 
on the BV formalism for Lagrangian systems [3], [1] with three major generalizations. The 
first generalization is that the consistency of gauge algebra is examined for the involutive 
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closure of the system of field equations, not for the action functional which might be even 
non-existent. The second is that the gauge identity generators are involved in the gauge 
algebra of the involutive closure independently from the gauge symmetries. The identity 
generators impose their own consistency conditions that are not identified by previously 
known method p] even in the Lagrangian case. The third is that the gauge algebra of the 
involutive system provides a convenient receipt ([8]), flHOj) for counting physical degrees of 
freedom. This formula is applied to the involutive equations in a covariant form, and it 
uniformly covers all conceivable instances. Let us also mention that the formula (IHOl) has 
been derived in the Appendix in a more general setting than it is actually utilized in the 
paper, as it also applies to the case of reducible gauge symmetries and gauge identities. 
Let us finally notice, that the involutive closure of field equations admits a BRST embed- 
ding along the lines of [25]. From the viewpoint of the corresponding local BRST complex 
|26] . the formula for the degree of freedom count (IE]), ( IHOl) can get a natural cohomological 
interpretation. This issue will be addressed elsewhere. 
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7. Appendix. 
Physical degree of freedom count 

In this Appendix, we explain the origin of formula (jH]) for counting of the physical degrees 
of freedom. 

The starting point for deriving the formula (IH]) is the notion of the strength of differential 
equations introduced by Einstein [H]. Roughly, the strength is a number that measures 
the size of the solution space. The "stronger" is the system of differential equations, the 
smaller is its solution space. It turns out that the numerical value of strength can be 
immediately related with the number of Cauchy data needed to define the general solution 
modulo gauge freedom, i.e., with the number of physical degrees of freedom. The original 
Einstein's argumentation was not mathematically rigor, its justification and explanation 
within the modern theory of formal integrability can be found in book [13]. For earlier 
discussions of the concept of strength of equation as well as numerous applications to the 
analysis of relativistic field equations we refer the reader to [I5], [IE], [IT], [T8] . 
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Let US first explain the Einstein's concept of the strength of field equations. Consider a 
set of fields 0* on (i-dimensional space-time with coordinates x'^. Assuming the fields to be 
analytical functions, we can expand them in Taylor series about some point Xq: 

oo ^ 

(58) <p\x) = J2 - ^orix - xor ■■■{x- xor^ . 

p=0 ^' 

Let Np denote the total number of terms of pth order in the expansion above. (The explicit 
expression for Np is given below.) As far as the fields 0* obey a system of PDEs, not all the 
Taylor coefficients V^^^.-.^p can remain arbitrary. Denote A^^ the number of monomials of 
order p that are left free in the general solution of the field equations. Obviously, Np < Np. 
On the other hand, if the field equations enjoy a gauge symmetry, not all the solutions are 
physically relevant; some of the monomials in the general solution come from Taylor series 
for the gauge parameters. Modding out by the gauge freedom, one can define the number 
A^' of gauge inequivalent monomials of order p entering the general solution. Now, the 
number of physical degrees of freedom per poin^ is given by 



(59) Af= lim 



p 



N'' 



p 



p^oo d — 1 Np 

This formula, that dates dates back to Einstein, defines the number of physical degrees 
of freedom as the growth of the number of "physical monomials" compared to the uncon- 
strained ones. 

Consider a system of PDEs 

(60) TM\d,f,...,d,,---d,^f) = 0, a = l,...,t, 

governing the dynamics of fields (f)\ i = 1,2...,/. The order of these equations equals to 
m. Substituting the expansion ( 158|) into the field equations (160|) and evaluating the result 
at X = xo, we get the system of algebraic equations 

(61) T,(y,Sy,;,...,y,;^...^,j = o. 

These equations follow from fIBU]) by simply replacing the partial derivatives of fields with 
the corresponding Taylor coefficients. In general, the solution space for these equations 
can be a very complicated algebraic variety containing strata of different dimensions. So, 
it might be problematic to choose the independent coefficients and compute their totality. 
But the task is considerably simplified if one considers the conditions on the higher order 
coefficients. After all, the lower order monomials do not matter for evaluating the limit 



'^Accordingly, the number of physical polarizations of the field is the half of Af. 
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( ]59ll . Differentiating ( 1601) k times by x's and setting x = xq, we obtain a set of algebraic 
equations of tlie form 

(62) J a^^ll■■■|lk ^^i-^m+k - = 0' 

where the functions J^'s and X's depend on y^'s of order less than m + k. Thus, for each 
given k, we get a system of linear inhomogeneous equations for the coefficients V^^j.-.i/^^j,- 
The matrix J' = J'k, defining the system, is called the symbol matrix of order k; it has the 
following structure: 



Here the round brackets mean symmetrization of the indices enclosed and the functions 
Ta{^) are given by the left hand side of equation (I6T1) . As is seen, the symbol matrix 
of order k is expressed in a very specific way trough the symbol matrix of order 0. The 
latter may be called the symbol matrix of the field equations f l60|) . For linear differential 
equations the symbol matrix Jq is just the highest-order or principal part of the system. 

We thus see that whatever the original system of field equations may be, there is an 
integer m such that the space of monomials of order p > m is determined by a finite system 
of linear inhomogeneous equations with coefficients depending on y^'s of order < p. The 
echelon form of the algebraic equations ( l62l) suggests to solve them one after another, so 
that at each step one deals with a finite system of linear inhomogeneous equations. This 
makes possible applying the usual theorems of linear algebra to evaluate the solution space. 

First of all, the number of linearly independent solutions to equations (!62|) crucially 
depends on the rank of the symbol matrix Ji. The symbol matrix in its turn is the function 
of the Taylor coefficients {v^^^...^^ }j=o constrained by the algebraic equations (EH), so that 
the rank of Jk can suddenly change. To avoid this complication we will restrict ourselves 
to those solutions of ( !6TI) for which the rank of the symbol matrix J'k{f) is maximal. 
This means that we consider only the general (opposite to singular) solutions to the field 
equations. Following [13], we will call {f^^j^...^j}Y=o principal coefficients, referring to 
the other Taylor coefficients as parametric. 

By Kronecker-Capelli's criterion, the system of linear inhomogeneous equations ( l62l) is 
compatible iff each left null-vector JC of the symbol matrix J'^ annihilates also the inho- 
mogeneous term X^. Clearly, the null- vectors of the symbol matrix, if any, can always be 
chosen to be functions of the principal coefficients alone. A crucial point is that the com- 
patibility criterion is automatically satisfied for differential equations in involution. The 
reason is very simple: vanishing of the function /C''^^ "^''Xa^i...^j. would otherwise give a 
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nontrivial constraint on the lower order coefficients that, in turn, would be manifestation 
of a hidden integrability condition. On the other hand, if the vector /C annihilates both 
the symbol matrix and the inhomogeneous term X^, then it defines an identity for the 
linear equations ( |62l) and this identity must follow from a gauge identity for the original 
differential equations ( 160|) . 

The general solution to an inhomogeneous linear system is given by its partial solution 
plus the general solution to the corresponding homogeneous system. In our case, the latter 
is completely determined by the symbol matrix. Thus, to evaluate the size of the solution 
space, we focus on the solutions to the homogeneous system. These form a linear space, 
whose dimension is given by the number of unknowns minus the rank of the symbol ma- 
trix. The number of unknowns {^^^^...y in (!62|) coincide with the number Np of lineally 
independent monomials of order p = m + k. It is easy to find that 

The rank of the symbol matrix jTfc can be computed as the difference between the number of 
equations (1621) and the number of left null- vectors of the matrix J7fc. The former is expressed 
through the binomial coefficients as 



(65) t 



k + d — 1\ ^fp — m + d—1 
I J \ p — m 



As was explained above all left null-vectors for the symbol matrix of involutive equations 
come from gauge identities. Each gauge identity 

(66) L'^Ta = 

is defined by differential operators 

n=0 

with coefficients depending on fields and their derivative up to some finite order j. If the 
highest coefficients ' are not all equal to zero identically, then the number q' is 

called the order of the gauge identity ( !T5|) . Differentiating ( 166|1 s times by x's and setting 
X = Xq, we find 

where dots stand for the terms involving <^'s of order less than m + s + q'. If s is large 
enough such that m + s + q' > j, then we are lead to conclude that the coefficients at 
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^9^^...^ + + / must be zero in order for the identity to hold. This imphes that the 

symbol matrix J's+q' admits the set of null- vectors {^Ai - a^} of the form 

Ai-As Ai As ' Ai-As •-^ai/ii---/i^_j^g/ 

All these null- vectors are linearly independent and their number is given by 

(68) (^ + ^-' 

In other words, each gauge identity of order q' for the involutive equations of order m gives 

(69) + /p-, + d-lA 

V / V p-Q J 

left null-vectors for the corresponding symbol matrix J'k provided that k is large enough. 
Here p = k + m and the number g = g' + m is called the total order of gauge identity. 

Now, suppose that the system (1^ involves equations of different orders: to equations 
of order zero (algebraic equations), ti equations of the first order and so on. Let us also 
assume that the system contains no hidden integrability conditions and becomes involutive 
upon adjoining trivial differential consequences. Then according to ( l65ll all these differential 
equations give rise to 



(70) X^t. 



n=0 



p — n + d — 1 
p — n 



linear equations for the parametric coefficients ^^^^...^^ with large p. Of course, only a finitely 
many terms are different from zero in the above sum. Let us further suppose that the field 
equations enjoy gauge identities of total orders n = 0, 1, 2, . . .. If all these identities are 
independent (irreducible), then the linear equations ( 16 2 p for the parametric coefficients of 
order p are possessed of exactly 

(71) ±4"-'"^"-^ 

dependencies (left null-vectors) provided that p is large enough. The difference between 
flTOj) and (I7T1) is the number of independent equations for unknowns V^^i-..^^- Subtracting 
this difference from we get the dimension of the solution space, that is, the number 
of independent monomials of order p: 



(72) A^; = / 



P J t-f, V p-^ 
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Now, we should take into account gauge freedom. Each gauge transformation has the 
form 

(73) 5,ct)' = Y,n'^--^-d,,---d,^e. 

n=0 

In this expression, the coefficients 7?.'s are functions of the fields and their derivatives up 
to some finite order and the infinitesimal gauge parameter e is an arbitrary function of x's. 
The number q is the order of the gauge transformation. The gauge invariance of the field 
equations (I6OII implies that 



(74) 5,T, = f/^n 

for some matrix differential operator U . Let us expand the gauge parameter in Taylor series 

00 ^ 

= - xqY' ■ ■ ■ (x - xoY- . 

n=0 

Then s-fold differentiation of the equality at Xq yields the identity 

{7K\ 7"A»l-"A'™.-T->iMm + l---A'm + 9c- _|_ . . . — fl 

{''->) ^ai fcAil---Mm + g + s "T — yi ■ 

Here the dots stand for e's of order less than m + g + s and all (/^'s are assumed to define a 
solution to the field equations. Since all e's are arbitrary, we conclude that the leading term 
in (175!) must vanish separately. This results in the set ^^^^,^^ "^1+^+=^ of right null-vectors for 
the symbol matrix JT^: 

These null-vectors span the space of dimension 

/m + s + q + d — 1\ /p + q + d — l\ 
\ m + s + q J \ p + q J ' 

In general, the system (160]) may enjoy several gauge symmetry transformations. We let r„ 



denote the number of the gauge transformations of order n. If all these gauge symmetries 
are independent (irreducible), then they make 



n=0 



p + n + d — 1 
q + n 



coefficients of {f^i - f^g} unphysical. Subtracting this number from ( !72|) . we get the number 
of "physically distinguishable" parametric coefficients of pth order. 
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Having computed Np and A''^ we are ready to evaluating the limit (l59ll . Making use of 
the asymptotic expansion for the binomial coefficients [15], [H] 

''p±n + d-l\ /'p + d-l\ ( n,, , ^/iM 

we find 



(77) N;:/Np={f-t + l-r)+^^^f2n{t^-L-r^) + o(\] . 
The numbers 

oo oo oo 

n=0 n=0 n=0 

coincide, respectively, with the total number of equations, gauge identities, and gauge 
symmetries. The leading term of the expansion 

(78) A = / - t + / - r 

is called the compatibility coefficient. Let us assume that the system (160|) is absolutely 
compatible, that means A = 0. Then comparing fl77|) with f lS^ . we finally arrive at the 
desired formula for the physical degrees of freedom 



(79) Af =^n{tn-l 



n ' n j 

n=0 



Vanishing of the compatibility coefficient A can be easily established under the assump- 
tion that each right null-vector of the symbol matrix J = J^^ originates from some gauge 
symmetry. To do this, we introduce the n x t-matrix 



Jaiip) — <Jai ^ Pm ' ' ' Pn 



m 5 



whose entries are polynomials in formal variables p^, fx = 1, . . . , d. It then follows from (|6 
that each gauge identity provides the left null- vector 

L^{p) = C''^--^^'p,,---p,^, 

for the polynomial matrix J{p), so that L°'{p)Jai{p) = 0. Similarly, each gauge symmetry 
transformation f lTS]) gives rise to the polynomial vector 

R\p) = n'^'-^''p,,---p,^ 

annihilating the matrix J{p) on the right, that is, Jai{p)R^{p) = 0. The vanishing condition 
for the compatibility coefficient flTHj) can be written as 
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It means that the rank of the rectangular matrix J{p), being computed by the number of 
left null-vectors, coincides with its rank defined in terms of right null-vectors. Clearly, this 
equality takes place for any matrix over an algebraic field, say M or C. It turns out that 
the same statement holds true for the matrices over the ring of polynomials in p's provided 
that the null-vectors L's and i?'s are linearly independent (over the ring of polynomials in 
p's) and span the right and left kernel spaces of the matrix J (see, e.g. [35]). 

There is also another, more direct, interpretation of the absolute compatibility condition. 
It can be shown [13j that the value f — t + I defines the number of arbitrary functions of 
d variables entering the general solution to the field equations ( l60l) . The equality A = 
/ — t + Z — r = then implies that all these functional parameters owe their existence to 
the gauge symmetries. To the best of our knowledge, example of field equations has been 
yet unknown that would not be absolutely compatible. Moreover, the results of [2] suggest 
that any system of ODEs is absolutely compatible and the same is true for two-dimensional 
field theories [36] • So, it is a very plausible hypotheses that any reasonable field theory is 
absolutely compatible. 

The above consideration can be extended to the field equations with reducible gauge 
symmetries and/or identities. Without further ado we just present the final formula for 
the physical degrees of freedom, which might be deduced by appropriate adjustment of the 
derivation in the irreducible case: 

(80) ^=11 ^(^n-(-ir(c+c))- 

m,n=0 

Here t„ is the number of the equations of order n; l"^ is the number of gauge identities of 
the total order n, and the reducibility order m; and is the number of gauge symmetry 
transformations of the total order n and reducibility order m. The total order of the 
generator of gauge symmetry/identity is defined inductively to be the sum of its order as 
a differential operator and the total order of a generator it annihilates. It is also assumed 
that the total order of the original gauge symmetry generators coincides with their order 
as differential operators, whereas the total order of gauge identities for the field equations 
is given by the order of the corresponding generators plus the order of equations they act 
on. For the field theories with irreducible gauge symmetries of the first order, equality (l80il 
was first derived in [18j . 

It is curious to note that the final formula ( IHOi) . being independent on (i, holds true for the 
one-dimensional systems as well, whereas the original definition (|59|) becomes meaningless. 
The proof of (IHOl) for d = 1 requires a different method, which is beyond the scope of this 
paper. 
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Let us concretize the formula ( |80l) for the special case of involutive Lagrangian second 
order equations as this case has a common interest in field theory. For the Lagrangian 
equations, the gauge symmetries and identities are generated by the same operators. The 
total order of the identities, however, is shifted by the order of the equations involved 
in, so the careful adjustment of the general relation (IHOl) for the second order involutive 
Lagrangian equations leads to the following count of physical degrees of freedom: 

(81) A/- = 2 ( t2 + Yl + 1)^1"^ 

\ n,m=0 

Here ^2 is a number of the Lagrangian equations, is the number of the gauge symmetry 
generators of the total order n and reducibility order m. In the irreducible case, (m = 0) this 
brings the well known relation for the degrees of freedom for the second order Lagrangian 
equations [19]: 

(82) Y = ^2 - 5^(n + l)r„. 

n=0 

J\f/2 has the meaning of number of physical polarizations. 
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